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We report a theory for analyzing nonlinear DC 
transport properties of mesoscopic or nanoscopic normal- 
superconducting (N-S) systems. Special attention was paid 
such that our theory satisfies gauge invariance. At the linear 
transport regime and the sub-gap region where the familiar 
scattering matrix theory has been developed, we provide con- 
firmation that our theory and the scattering matrix theory 
are equivalent. At the nonlinear regime, however, our theory 
allows the investigation of a number of important problems: 
for N-S hybrid systems we have derived the general nonlin- 
ear current-voltage characteristics in terms of the scattering 
Green's function, the second order nonlinear conductance at 
the weakly nonlinear regime, and nonequilibrium charge pile- 
up in the device which defines the electrochemical capacitance 
coefficients. 
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I. INTRODUCTION 

The physics associated with quantum conduction in 
various low dimensional hybrid superconducting systems 
has been a major focus of research at presentiJa. Due to 
advances in controlled crystal growth and lithographic 
techniques, it is now possible to fabricate various submi- 
cron hybrid superconducting structureSpWhere accurate 
experimental measurements can be madeO. An important 
theoretical task is to be able to predict, for general meso- 
scopic or nanoscopic hybrid systems, transport properties 
such as the nonlinear current-voltage characteristics and 
the nonequilibrium charge distribution inside the system 
as a function of the applied bias voltage. Our theoretical 
understanding of quantum transport in these very small 
N-S hybridj-sijstems has been achieved by scattering ma- 
trix theorya'cl a|nd_.by non-equilibrium Green's function 
(NEGF) theoryO^y. 

To analyze nonlinear transport coefficients, i.e. coeffi- 
cients which appear in front of nonlinear powers of bias 
voltage, in principle one must make sure gauge invariance 
of the theory. This means that theoretical results should 
not change when bias voltage applied at all the device 
leads is changed by the same amount. This is a nec- 
essary condition for any tijaiisport theory and has been 
recognized in the literatureEJ. Consider a device which is 
connected to the outside world by several leads a where 
bias voltage Vq is applied. When Va ^ Va + v where 
w is a constant, the calculated results (such as current) 
will not change if the electrostatic potential U inside the 



device is also changed by the same amount v. However 
U = U{r) which is the Hartree potential, can only be 
obtained by solving a self-consistent problem. In other 
words, to satisfy gauge invariance one is necessarily re- 
quired to consider Coulomb interactions at least at the 
Hartree level. Furthermore, in general when external bias 
voltage is applied to a device, the flow of charge carriers 
through the device could polarize the system due to long 
range Coulomb interactions. For a macroscopic metallic 
conductor, the polarization can be safely neglected since 
interaction is well screened. However for mesoscopic scale 
and nano-scale conductors the polarization could be very 
important. This also requires self-consistent analysis. 

For I nor mal conductors, Biittiker and his co- 
workergl2lll3 developed an approach based on the scat- 
tering matrix theory to deal with the second order non- 
linear conductance coefficients. This theory can also rbe 
extended to higher nonlinear order in DC situationsli3. 
On the other hand for a N-S hybrid system, despite the 
many theoretical investigations on its quantum transport 
propertyO, the important issue of gauge invariance has 
not been clearly addressed so far. In light of this un- 
satisfactory situation, in this paper we report the de- 
velopment of a proper nonlinear transport theory which 
satisfies gauge invariance for mesoscopic or nanoscopic 
N-S hybrid device systems. Our theoretical formulation 
is based on nonequilibrium Green's function approach 
where the quantum transport problem is solved in a self- 
consistent manner. We have derived analytical expres- 
sions for the general I-V characteristics including the 
sub-gap behavior of the N-S device in terms of Green's 
functions which are numerically calculable. In the weakly 
nonlinear regime we further derived the second order non- 
linear coefficient by solving the characteristic potential. 
Since our theory can deal. with charge polarization in the 
presence of transportlijO, we have also derived the lin- 
ear and second order nonlinear charge distributions which 
define the electrochemical capacitance of the system. Fi- 
nally, at linear order for which scattering matrix theory 
(SMT) has been developed for N-S systems, we prove by 
explicit calculation that our theory and the SMT theory 
are equivalent. This way the widely used SMT theory 
and the NEGF theory are formally connected for N-S 
systems through a Fisher-Lee relation which we derive in 
this work. 

The rest of the paper is organized as follows. Section 
II and III present the gauge invariant nonlinear trans- 
port theory for the N-S system and confirming its equiv- 
alence to the familiar scattering matrix theory at the 
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linear transport regime. Section IV provides several ap- 
plications of our theory. Section IV is a short summary 
for this work. 



II. GAUGE INVARIANT THEORY 

In this section we develop the gauge invariant nonlinear 
transport theory based on NEGF for N-S hybrid device 
systems. To be specific, the N-S system we consider is 
a quantum well connected to a normal metal lead and 
a superconductin£- lead. It is described by the following 
HamiltonianQO'Ll in the second quantized form, 



where 



H — Hl + Hri + Hd + Ht 



(1) 
(2) 



k.a.i 



(3) 



(4) 



+ tR 



-ievRT J 



(5) 



Here Hl describes noninteracting electrons in the left 
lead (assumed to be a normal- metal), a]^ ^.^ {o-L.kcr) is the 
creation (annihilation) operator of electrons in the left 
lead, and vl is the bias voltage applied on the left lead. 
Hr describes the right lead which is a superconducting 
lead whose gap energy is A. Hd is the Hamiltonian of the 
quantum well with multiple discrete energy levels char- 
acterized by index i and spin a. Ht denotes the coupling 
part of the Hamiltonian, ta^ki {oi = L, R) is the hopping 
matrix and for simplicity is assumed to be independent 
of the spin index a. In order to obtain the Hamiltonian 
(^ and (H) , we have performed a unitary transformation 
similar to that described in Rcf. ^, so that the bias volt- 
age of the right lead, vr, appears as a phase factor in the 
hopping elements. Finally, it is important to note that 
we must include the internal Coulomb potential U inside 
the scattering region for further analysis (see below). 

Our analysis follows that of Ref. ^ by applying the 
NEGF theory, iterating the equation of motion for the re- 
tarded Green's function, and applying the Keldysh equa- 
tion for the lesser Green's function at equilibrium. The 



current flowing through the normal lead is derived to beE 

{e = n = i), 



I = lA + h 



(6) 



with 



Ia= 2 j ^Tr [TlGI^TlGI^] [fL{E + vl- vr) 
-hiE-VL + VR)] , (7) 

— PK(i?)Tr [TlGI^TrGI, + TlGI^TrG'^^ 



-J^^(^lGIiTrGi2 + ^lG12^rGii) 
MfLiE + VL-VR)~ fR{E)] 



(8) 



where Gn and 6*12 are the matrix elements of the 2x2 
Nambu representation.. Here an important departure 
from previous theoryQ is the explicit inclusion of the 
Coiulomb potential t/(r) in the Green's functions G^i and 

Gl,{E) ^[E-Hd + U-VR--El, 

- S^^A-S^^i]-! (9) 



^[E + h:i+vr-U-^'22] ■ (10) 

Once the electron and hole Green's functions G\i and A'' 
were obtained, GI2 is calculated by 

G[2 = Gl.Sl^A^ . (11) 

The self energy S*" = S£ + SJj, is derivedi to be 



^1 = 











•pa 



(12) 



where EJ^ = Pq, — iTa/2 is the self energy of the lead a 
in the normal case. Here Pa is the real part and Fq-Is 
the linewidth function. In the wide bandwidth limitEj, 
the self energy for the superconducting lead is 



2 1/32 



P2 



PI 



(13) 



where (3i = ^^/^/E^^ - A^, P2 = A/^/E^ - A^, n = E 
when < A and k = otherwise so that the electron- 
hole symmetry is preserved. The dimensionless BCS den- 
sity of states (DOS) is given by pR = e{\E\ - A)/3i. 

We emphasize that the crucial step in developing the 
gauge invariant nonlinear DC theory is to include the 
internal potential lanjdscape C/(r) into the Green's func- 
tions self-consistentlyll3. In this work we deal with it at 
the Hartree level, hence t/(r) is determined by the self- 
consistent Poisson equation 
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V^U{x) = ^47Ti{G<iE,U)). 



(14) 



where Gf^ is the electron lesser Green's function in real 
space and x labels the three dimensional position. From 
Ref. ^ Gil determined by iterating the Keldysh equa- 
tion at equilibrium, the result isB 

— [GliTLGtifLiE + VL-VR) 
+ Gl^TLG^^jLiE-VL+VR)] 

/dE 
—pj,{E)fR{E) [GIiTrGIi + GI^TrG^^i 



+ I 



(15) 



Eqs. (0), (H), and (|^ completely determines the non- 
linear I-V characteristics of the N-S hybrid system: they 
form the basic equations of the gauge invariant nonlin- 
ear theory. The self-consistent nature of the problem is 
clear: one must solve the quantum scattering problem 
(the Green's functions) in conjunction with the Poisson 
equation. It is easy to prove that the current expression 
Eqs.(^ and (||) are gauge invariant: shifting the bias 
potential everywhere by a constant w, — > -I- w so 
that U ^ U + V (since the boundary condition of U is 
changed by w), / from Eqs.(|7|) and (^) remains the same. 
Eqs. (|^), (^), and ( p^ also form the basis for numerical 
analysis of I-V curves for the N-S system. For instance 
one can compute the various Green's functionS|Xj and the 
coupling matrix F using tight-binding modelal, and the 
Poisson equation can be efficicntljj solved in real space 
by powerful numerical techniquealj. 

In the simplest approximation, the gauge invariant 
condition can be satisfied by putting a gate voltage Vg 
as was done in Ref. ^ so that one treats the system as 
a three probe conductor with external voltages Vl, Vr 
and Vg applied at the probes. In general, the internal 
potential is a nonlinear function of Va (see section IV for 
details), but as a first approximation one expands it in 
terms of Va in the small voltage limit. 



U = UlVl + UrVr + UgVg 



(16) 



where Ua (r) is the characteristic potential which satisfies 
the sum rule Ua = 1. If one makeaJ a further approx- 
imation by assuming ul = ur = 0, the sum rule gives 
Ug = 1 and U = Vg, i.e. t/ is just a constant potential 
shift under these approximations. 

In distinct contrast to the constant U model, the the- 
ory presented in this section is a microscopic gauge invari- 
ant theory. Furthermore, in order to discuss charge polar- 
ization and electrochemical capacitance in the presence of 
transport, one has to include the self-consistent Hartree 
field rather than just include a constant gate voltage: 
one can easily confirm that the constant U model corre- 
sponds to the local charge neutrality approximation, it 
will therefore not give rise to any charge polarization. 



III. THE SCATTERING MATRIX AND 
FISHER-LEE RELATION 

The Fisher-Lee relation!^ relates scattering matrix of a 
conductor to its Green's function. This relation has been 
widely used in analyzing quaatum transport through 
mesoscopic normal conductorsEll. To make a formal con- 
nection of our theory presented in the last section and the 
scattering matrix theory for N-S systems, it is necessary 
to derive a Fisher-Lee relation for the N-S system. In this 
section we provide such a derivation which allows us to 
make comparison between our theory and SMT. Because 
results (gauge invariant) of SMT for N-S system are only 
known to linear order in bias, we reduce our thecay to 
that limit in this section. In addition, the SMTE3i3 has 
so far dealt with the sub-gap region {\E\ < A) hence we 
will focus on the Andreev current Ia only, although our 
theory is applicable beyond this limitation as shown by 

Eq. &. r-. 

To prove that our theory reduces to that of the SMTE3, 
we go to the wide bandwidthcj Jimit at zero bias within 
the gap. In this limit we havelia = —iTa/'2, hence 
Eqs. (|, |l^, 0) become 



Gil 



Gr 
22 — 



E + H*+iTL/2 + irRPi/2 
X 



E-Hd + 2TL/2 + iTRl3i/2 
X 



and 



where 



12 — '-^21 — 



iTRfi2/2 
X 



X={E + H*+ irL/2 + tTRj3i/2){E 
+ iTL/2 + tTR(3i/2) + l3lTl/4 . 



(17) 



(18) 



(19) 



(20) 



Note that these reduced Green's functions (|T^,|l^,^9|) can 
be collectively defined by 
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E - Heff +irL/2 + irRpi/2 
for which the effective Hamiltonianii is given by 



Hd 

iTrI32/2 



iTRP2/2 
-H* 



(21) 



(22) 



The scattering matrix approach of Ref. ^ is recov- 
ered by letting E ^ (which is equivalent to setting 
the parameter of scattering theory a — —i, see Ref. p2l ), 
therefore /?i = and (32 = — «• We thus obtain 



Hd rR/2 
r/?/2 -H* 



(23) 
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which is exactly the same as that of the scattering ma- 
trix theory of Ref. We thus conclude that the wide 
bandwidth limit of our theory is equivalent to SMT at 
zero bias. 

The equivalence of our theory and SMT can be further 
discussed throug|ti-|the Fisher-Lee relation. This relation 
was first derivedcj for normal conductors, and here we 
derive it for N-S systems within the sub-gap region. For 
simplicity we assume that the scattering matrix in the 
normal side of the N-S system is of the Breit-Wigner 
form: 



Safi = 



E-Ea+iTl2 



(24) 



with F = Fi + From Ref. g, the scattering matrix of 
the N-S hybrid system is given by 



with 



and 



She{E) = ae-"l'sl^{-E)M,S2i{E) 
M, = [l-a\s22{E)sl2{~E)]-^ 
E - VE^ - A2 



Using Eq.(p4|) we obtain 

S22{E)S*22{-E) = 1- 



-F1F2 + 2iET2 



iE-Eo + iT/2){E + Eq+ iT/2) 



and 

S2l{E)sl2{~E) 

Hence 



riF2 



{E-E0+ iT/2) {E + E0 + iT/2) 



She = ae 



Y 



where 



Y = (1 - a^){E -Eo + iT/2){E + Eq + iT/2) 
- a^{T{T2-2iET2) . 



Finally, note that I - = 2ay/E'^ - A^/A, we ob- 
tain a relationship between the scattering matrix She and 
Green's function 

sUE) = ^^^0t=e-^%E -Eo + i^)iE + Eo + ^^) 



Tl/A+i2E + iTi){l + a 



y/E^ ~ A2 ' 2 



riF2A 

2y/E^ - A2 
rLe-'*Gl2 



e-"^/X 



(25) 



where we have used the fact that 1 -I- a A / V-E^ — A^ = 
i;/V^2_ A2, Ti = Fi, r2 = Ffl, and-Eq. ^ for 
G\2. Eq.dl) is the Fisher-Lee relationEa for the N-S 
system. Using Eq. ( p5| ) it is straightforward to show that 
when the Fermi energy is inside the gap, Eq.(^ gives 
the same result as that of SMT in Ref. ||J| with a = 
exp(— j arccos E/ /S). 



IV. APPLICATIONS 

In this section we present detailed analysis for a num- 
ber of situations where analytical expressions can be ob- 
tained in closed form. These are resonance Andreev re- 
flection coefficient at the linear regime; the second order 
weakly nonlinear conductance; and the nonequilibrium 
charge distribution. 



A. Resonant Andreev reflection 

The phenomenon of resonant Andreev reflection^! is 
the situation where the Andreev current is dominated 
by a resonance transmission through a level Eq inside the 
quantum well. For this case the Andreev reflection coeffi- 
cient Ta{E) = rr[FLG^2rLG"2] can be derived assuming 
a Breit-Wigner form of the scattering matrix (|2j). From 
Eq. ( ^5|) and note that at sub-gap region energies Eo and 
E are small therefore we can take a ~ —i, we obtain 



ir2 



2{E^ - El) - (r2 + JF2)/4 + iE{T + 5T) 



where (5F = Fi — r2. After simple algebra the Andreev 
reflection Ta is foimd to be 



Ta 



J- 1^ 2 



4(^2 -El + r(5F/4)2 + TlTl + £;2(r2 + ^r2) 



(26) 



This is the Breit-Wigner formula for N-S system. Note 
that this expression is different from the formula of Ref. || 
which is only valid at = 0. Several interesting observa- 
tions of ( p6| ) are in order. First, let's consider resonance 
energy being at Eo = 0. For this situation we have 



Ta 



p2p2 
^1^2 



4(i;2-HF(5F/4)2 



TlTV 



Therefore: (i) if Fi < F2 or ^F < 0, there are two peaks 
with Ta ^ I a.t E ^ ±V-F5F/2; (n) when Fi = F2, we 
have Ta = TlTl/[AE'^ + FfF|], so that there is only one 
peak with Ta = I aX E = 0] (iii). for Fi > F2, there 
is only one peak with Ta = F?F|/[F2(5F2/4 + F?F^] < 1 
a,t E = 0. Hence when energy level of the quantum 
well is aligned at the center of the superconducting gap, 
the resonance Andreev reflection is characterized by one 
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or two peaks depending on how the normal and super- 
conducting leads are coupled to the quantum well. Sec- 
ond, when Eo is nonzero, we have: (i) Fi < two 
peaks with Ta = TItI/[TItI + El{T'^ + SV^)] to occur at 
=Ei~ T6T/4:. (ii). Ti > Ta, just one peak at E = 
if r5r/4 > E^, otherwise two peaks at E^ = E^-TST/A. 
Note that when Eo is nonzero, the peak values are always 
less than one. 



B. Weakly nonlinear regime 

For weak nonlinearity we can expand all quantities in 
terms of the external bias voltagetll and obtain results 
order by order. Such an expansion makes sense when bias 
is finite but small. TJais approach was adapted in SMTtll 
and response theorytj for analyzing normal mesoscopic 
conductors. For the N-S system we will derive formula 
for the local density of states (LDOS) and the second 
order weakly nonlinear DC conductance. These are the 
interesting quantities for weakly nonlinear regime. 

In both SMllB and response thcorytj, LDOS plays a 
very important role. From our NEGF theory LDOS can 
be easily derived from the right hand side of Eq. (14), 
which is the charge density, with the help of Eq. (15). 
Here jwe shall present the explicit expression at the lowest 
orderfij expansion in external bias. Hence we seek the 
solution of U (r) in the following form, 



a Q/3 



(27) 



Pinj = {due/dE - dnh/dE){vL - vr) 
- {l/2){d'^nJdE^ + d^nh/dE'^){vL 



VR? (30) 



is the injected charge from the normal lead. due/dE is 
the injectivity of electron, i.e. the DOS for an electron 
coming from left lead and exiting the system as an elec- 
tron, 



dn,{x)/dE - j {dE/27r){-dEfL){G[,^o^LG'i,^o)^ 



(31) 



In addition drih/dE is the injectivity of a hole, i.e. the 
DOS for a hole coming from left lead and exiting the 
system as an electron, 

dnh{x)/dE - [ {dE/2n){-dEfL){G\^,aTLGl^,a)^., . 



(32) 

Finally eP'njdE^ is the derivative of dn/dE with respect 
to energy. Note that Eqs.(^l]) and ( p2[ ) are the same as 
that defined ia the scattering approach of Gramespacher 
and BiittikeiEl. 

In Eq.(|29|) the induced charge due to long range 
Coulomb interactions is derived to be given by 

p,Ux) = -J {dE/TT)fLlm [G^i,o(«i - SJa^^L^SS^i) 



(33) 



where (.r) . and Wa/3..(r) are the characteristic 
potentialsEi[E3. It can be shown tiLat-the characteristic 
potential satisfies many sum ruleslliliij, Ua = 1 and 



0, where the subscript {/?}; is a short no- 



tation of / indices 7|r^, i],- ■ ■. Expanding Gfi of Eq. (|l4| ) 
in power series of WqE3, we can derive equations for all the 
characteristic potentials. In particular the expansions are 
facilitated by iterating the following Dyson equation to 
the appropriate order 



A'' = AS - {vR 



U) Al 



and 



11 — "-^11,0 ~ "-^11 \^ 
^21 ^11,0 



VR) G^i,o + Gn (A'- - Al) 

(28) 



where Aq and G\iq are equilibrium hole and electron 
Green's functions. The expansion of G\2 can be made 
similarly. At the lowest order, we thus obtain the local 
charge density in the presence of transported. 



where H is the general- 

ized Lindhard function which-pdjices to the Lindhard 
function of normal conductortZHiillij in the limit A ^ 0. 
For example, using the wide bandwidth limit expressions 
(such as Eq.(p7|) for the Green's function, or using the 
Breit-Wigner form Eq.(p4[) for the scattering matrix, at 
small Fermi energy (so that (3i — Q and [32 = —i) we 
can calculate the Lindhard function H exactly at zero 
temperature from its definition above: 



H 



Li R 



TT 2^2 

arctan 

2 



(r 



r|j)/2' 



Li R 



(34) 



where we have set the quantum well level Eq = Q. Hence 
the Lindhard function is a smooth function increasing 
with energy monotonically. 

With these quantities the Poisson equation becomes, 

dne{x) dnh{x) 



dE 



dE 



(35) 



p{x) — i{Gfi — GflQ)xx — Pinj + Pi; 



where 



(29) 



'Ull{x') = 47r( 



dhf.(x) dnh{x) 



dE 



dE 



(36) 
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where dheJdE, and dfih/dE are the second order 
injectivitiestSH. These partial differential equations can 
at least be solved numerically. However to avoid nuniMi^. 
ics one may apply the quasi-neutrality approximatior£3 
by neglecting the spatial derivative in Eqjmq) . This way 
the characteristic potential is obtained aa^J 



UL 



Arip duh , , 

— -) n 

^ dE dE " 



(37) 



In terms of the characteristic potential we now derive 
the second order nonlinear conductance due to Andreev 
reflection. In the weakly nonlinear regime, only the An- 
dreev current I a is relevant which can be expanded in 
terms of external bias voltage difference v = vl — vr, 

Ia — Giiv + Gixiv^ + ... 

From this definition of conductance coefficients Gn and 
Gill, we expand Eq.(|^) in terms of v to obtain, 

Gii^^j {dE/2Ti){-dEfL)TA 



and 



Gil 



dE 
2^* 



'dEfL)Tr 



dGA 
dU 



■UL 



(38) 



where Ga = FlGijElGij and dGA/dU is easily calcula- 



ble using Eq. (|2^) and the relation in Ref. To compare 
with the second order conductance of normal conductor 
Gf^i, we note that Gf[i has two contributionsE£lE3. One 
of them comes from Coulomb interaction, 



Gfii = J {dE/2TT)Tr[G^{rLGl,UL + ULG^ri 
l/2rLGS - l/2G^,rL)G'oTj,]dEf 



(39) 



However, for NS system, if the Coulomb interaction is 
not important (when ul = 0), we would have Gm = 0. 
For example, for an ideal ballistic wire, or for a symmet- 
ric quantum well at a resonant tunneling point, every 
incident charge is perfectly Andreev reflected. There- 
fore, for these examples no charge accumulation is pos- 
sible. From Eqs.(|l]) and (p^), we can easily verify that 
drie/dE = dufi/dE near a resonant point for a symmet- 
ric system and hence a vanishing Gm since ul = from 
Eq.(37). In contrast, when = 0, Gf^i is nonzero from 
Eq.(39). 



C. Electrochemical capacitance 



Qi = Ciiv + -jCwiv^ + •■■ = C{v)v 

this defines the electrochemical capacitance coefficients 
Gil, Gill and the general voltage dependent electro- 
chemical capacitanceE2l C{v). It is not difficult to confirm 
that the first two coefficients are 



Gil = 



dUe dUh 

"^""^lE - -dE ^ 
dxdx'Il{x, x')ul{x') 



(40) 
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dx{ 



drip 



dfih 



dE dE 



) 



— f dxdx'Il{x,x')uLL{x') 



(41) 



To get some physical insight for these coefficients we 
consider the discrete potential modellij. In addition, we 
parameterize the characteristic potentials by the geomet- 
ric capacitance Gq, in terms of which the Poisson equa- 
tion becomes 

Coiui - uii) = {Dj - D'})v - Hiui = Ciiv (42) 

- C„{uj ~ uji) = {Djj - D'}j)v - Uiiuii (43) 

where we have set D"^ = Jj dx{dne/dE), = 
J J dx{dnh / dE) , H/ = dxli{x, x) and small bias limit is 
assumed. We solve the characteristic potentials ul and 
ull through these two equations in terms of Gq. This 
leads to the following expression for the electrochemical 
capacitance coefficient Gn for a N-S system: 



Gil = 



Co' 



h; 



H 



(44) 



In particular, in the limit of gap A ^ 0, from Eqs.(^l]), 
and dl) we obtain = 0, U ^ dn/dE, and 
— driL/dE where dn^/dE is the injectivity of left 



D 

lead. In this situation Eq. (|44|) reduces to the expres- 
sion of tha electrochemical capacitance for a normal 
conductoiEj. Let's consider a symmetric tunneling N- 
S system. At the resonant point, the electron will be 
reflected as the hole due to the Andreev reflection. As a 
result, the capacitance Gn vanishes since D"" = D'^' and 
there is no charge accumulation. 



Using the NEGF theory one can also investigate the 
nonequilibrium charge distribution inside the N-S sys- 
tem. For this purpose we divide the system into two 
regions: in region I the charge is positive and in region II 
it is negative. The total charge in region I can be calcu- 
lated using Eq.(^9|): Qj = Jj p{x)dx. Expanding Qj in 
powers of v in the following formEj£3 



V. SUMMARY 

In this work we have developed a gauge invariant 
NEGF theory for hybrid N-S systems. This theory ex- 
plicitly takes into account the long range Coulomb in- 
teraction in the normal region. When the Fermi energy 
is inside the gap, we have shown that at linear regime 
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this theory is equivalent to the scattering matrix the- 
ory. We have also derived the Fisher-Lee relation for 
the N-S system hence a formal connection between our 
NEGF theory and the SMT is made for these systems. 
Because of gauge invariance, our theory is applicable for 
nonlinear regime for which we have derived an explicit ex- 
pression for nonlinear current-voltage characteristics for 
N-S devices. This result can be further simplified in the 
weakly nonlinear regime, for which we have analyzed the 
second order nonlinear conductance and the generalized 
Lindhard function. It is interesting the see that for N-S 
systems the concept of injectivity is naturally extended 
to include the injectivity of holes: these quantities auto- 
matically appear in our formalism. Our theory included 
charge polarization effect hence can be applied to analyze 
the electrochemical capacitance coefficients at the linear 
and nonlinear orders in bias. In particular we have de- 
rived an analytic expression of the linear electrochemical 
capacitance of N-S system within the discrete potential 
model. 

While this paper concentrated on the development of 
a theoretical formalism in terms of the Green's func- 
tions, it is obvious that numerical computations can be 
carried out applying the analytical expressions derived 
here. This way one can avoid the various approxima- 
tions used here in order to obtain closed form results. Of 
particular interests are the investigation of nonlinear I-V 
curves without the wide bandwidth approximation; the 
calculation of nonlinear conductance coefficients with- 
out the quasi-neutrality approximation; and the study 
of nonequilibrium charge distribution without the dis- 
crete potential approximation. These, however, will be 
the subject of a future report. 
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